Device-independent entanglement quantification 
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We present a method for quantifying both bipartite and multipartite entanglement in a device- 
independent manner. Specifically, from the observed correlations or the amount of violation of 
a given Bell inequality, we present a lower bound on the amount of entanglement present in the 
underlying system, as measured by the negativity or the genuine negativity. Some of these bounds, 
such as the one we obtain for the Clauser-Horne-Shimony-Holt Bell inequality, are shown to be tight. 
These bounds provide further information about the dimension of the underlying Hilbert space or 
the type of entanglement responsible for the observed correlations. Our method relies on a novel 
technique to characterize correlations derived from quantum states having additional properties, 
such as a positive partial transpose. 

PACS numbers: 03.67.Mn, 03.65.Ud, 03.67.a 



Introduction. — The analysis of most quantum experi- 
ments relies on a preeise eharaeterization of the equi 
ment employed to generate the data, e.g., see Ref. 
and references therein. However and maybe quite sur- 
prising at first, many tasks in quantum information pro- 
cessing can be achieved without sueh eharaeterization or 
knowledge, allowing the conclusions drawn to be device- 
independent. This approach first started in quantum key 
distribution where security statements that are indepen- 
dent of the precise implementations are clearly desirable. 
Indeed, already the security of the Ekert protocol Q 
relies on a Bell inequality violation Q rather than on 
non-distinguishability of special quantum signals. This 
direction was further motivated by device testing scenar- 
ios where both the true device and the test apparatus 
come from an untrusted source 0, Hf. In recent years, 
a variety of other device- independent quantum informa- 
tion tasks has been proposed, including further ciuan- 
tum key distribution [6|-l9| or self-testing schem es l id . Ill 



tests for different types of entanglement [12 



...... ...„_^..^^ ^ ^j, ways 

to certify the state space dimension |17H19| . the nature 
of measurements (2^ 21 1 as well as randomness expan- 
sion [g^]. Some of these examples are only partially 
device-independent but follow the general spirit of re- 
moving assumptions about the quantum devices. 

In this paper we introduce a technique to characterize 
the correlations obtained by measuring quantum states 
which are positive under partial transposition (PPT) [2^ . 
This sets a framework for device-independent entangle- 
ment quantification in the bi- and multipartite case. 
Moreover, this technique provides information on both 
the dimension and the type of entanglement present in 
the measured system. Note that such a quantification 
has so far only been achieved for the simplest experi- 
mental scenario together with a qubit assumption |24{ or 
even further measurement properties 25|. We employ 
the negativity (26| as our primary entanglement mea- 
sure [27[ . This provides further insights into the bipartite 



Peres conjecture [281, which implies that no PPT entan- 
gled state can violate a Bell inequality. Here, we show 
that a violation by any such state can at best be very 
small, if not vanishing, for a large class of Bell inequali- 
ties. 

Problem definition. — Let us start by considering a bi- 
partite Bell-type experiment where each party can em- 
ploy different measurement settings x, y giving rise to 
outcomes a, b, respectively; we denote the resulting con- 
ditional probability by P{a, b\x, y). This distribution has 
a quantum representation if there exists a quantum state 
Pab and local measurement operators Ma\x and M^y 
such that P{a,b\x,y) = tr(p^BMaj^ (g) M^Ij,). Our aim 
is to characterize the set of P{a,b\x,y) obtainable from 
Pab satisfying certain properties. If pab is only required 
to be a valid quantum state, we recover the original ques- 
tion leading to Tsirelson's bounds [29l - f3ll |. But one can 
require pab to satisfy more constraints, such as being 
PPT, i.e., that its partial transpose 23 1 p^^ has only 
non-negative eigenvalues, p^g > 0. One can also ask 
for a characterization of P{a, b\x, y) if the amount of 
entanglement is bounded E{pab) < -E'max using a pre- 
ferred entanglement measure E. This characterization 
task generalizes naturally to the multipartite case, e.g., 
to describe if the tripartite distribution P(a, 6, c\x, y, z) is 
quantum, biseparable (l^ . [sH or originates from a PPT 
mixture [33|. 

Our solution is a superset characterization, analogous 
to the method of Navascues-Pironio-Acin (NPA) [s^ : For 
instance, in the bipartite case we show that a given dis- 
tribution P — P{a, b\x, y) can only originate from a PPT 
state if a special matrix x[P-i ''At that linearly depends on 
P and on some unknowns u, satisfies x[P, m] > and 
x[P, u]^"* > 0. If it is impossible to find such parameters 
M, then P has no PPT quantum representation. 

Note that the aimed characterization must hold for all 
possible dimensions, measurements and states with the 
desired property. However, since any measurement op- 
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erator corresponds to a projector in higher dimensions, 
we can assume the projection property in the foUowing, 
i.e., the relation Ma\xMa.i\^ = (5aa'-^a|a: for the operators 
M(j|^ on system A for all a;, a and a' . This follows from 
Naimark's extension [sJl that also preserves any entan- 
glement monotone. Also, we shall simultaneously employ 
the notations M^^^^ and Ai for measurements operators 
on system A, likewise for other systems. The set {At} 
contains the identity operator Aq ^ I and all but one 
measurement operator Ma\x for each setting. Hence one 
has the aforementioned projection property and an iden- 
tity relation AiA^ = AqAi = Ai for all i. 

Technique. — In order to characterize correlations in re- 
lation with entanglement we combine the concept of ma- 
trix of moments 35 - 3^ with the special partial informa- 
tion of the device- independent setting [sO, [sH 113, IHj . 

Let us start with the matrix of moments for the bipar- 
tite case and consider first the scenario where the state 
Pab and measurement operators Ma\xT J^b\y ars known. 
To this scenario we associate two completely positive 
(CP) local maps Aa, As that we apply to the quantum 
state x[p\ = x[pab]ab = Aa «) ^^b[pab]- Here A and 
B denote the respective output spaces. Specifically, con- 
sider the local map Aa [p] = X) P Kn where the Kraus 
operators arc given by Kn = J2i K)AA("Mi' I^)a' 
\i) A arc orthogonal basis states of Ha and Ha respec- 
tively. Using a similar map for B one obtains 



x[p] = 



ijkl 



\ij)ab {kl\tY[pABAlA, 



b}b,]. (1) 



Thus the matrix x[p\ is just a matrix of certain expec- 
tation values. Since the local maps can also be defined 
using higher moments, e.g., by choosing Kraus operators 



. Ail, , we shall 



refer to % as a moment matrix of level i if it contains all 
£-fold products of Ai. Since both sets {Ai}, {Bj} con- 
tain the identity, the trace of the underlying state is a 
matrix entry that we refer to as x[p]tr = tr[p]. Finally, 
note that by the structure of these local maps we have a 
couple of important relations: i) if p > then x[p] ^ 0, 
ii) if p^-* > then x[p]^'* > 0, and iii) if p^-'' > then 
xIp"^"*] ^ 0. This matrix of moment approach can analo- 
gously be defined in the multipartite case. 

A device-independent characterization draws conclu- 
sion only from the observed correlations, hence, many of 
the entries of x are unknown a priori. However, even 
without this information the matrix x[p] has a structure 
which follows from known relations that hold indepen- 
dently of state and measurements: 1) Ai,Bj are Hermi- 
tian operators. 2) Ai,Bj satisfies the above mentioned 
projection property and the identity relation, 3) cer- 
tain entries correspond to the observations P{a, b\x, y) = 

tv[pABM^yx®M^y\. 

Via this partial information we can decompose without 



loss of generality each matrix of moments x[p] as 

x[p]^x[PM^x''^{P) + x°''''\u) 

= ^ P{a,b\x,y)Fabxy + ^UyFy, 



(2) 



a,b,x,y 



i.e., into one fixed part that linearly depends on the ob- 
served data x^^{P) — '^P{0',b\x,y)Fabxy and into an 
orthogonal, open part x°^''"(w) = Su'^f-^'i' which would 
be known only by the knowledge of state and measure- 
ments. Here all operators = are Hermitian. Note 
that the constraint x[p]tr = xi-Pj^'ltr = 1 is fulfilled auto- 
matically if the probabilities P are normalized. We give 
an example how the relations 1) — 3) provide the form 
given by Eq. ^ in the appendix. 

Connection with the NPA hierarchy. — At this point 
we like to connect the present technique to that of 
NPA [s^, 41 1, the best known method to characterize 



quantum correlations. For their method, one can identify 
a likewise construction x^^^[p] = A[pa_b], but with A be- 
ing a global CP map which already ensures that if p > 
then x^^^[p] ^0- K one uses the operator-sum ansatz 
X^^^[p] = Em I^mpL]n where = J2s !«) ('"I Os with 
|m) , \s) being respective basis states for the bipartite 
in- and output Hilbert spaces, this leads to x'^^^W = 
J2 \s) {t\ tr[pABOlOs]. If this operator set {Os} consists 
of all £-fold products of measurement operators, then im- 
posing the constraint x'^^^W ^ corresponds to the 
£NPA_^^ step in their hierarchy. 

Therefore a bipartite moment matrix x of level £ as de- 
fined above and a 2^-step x^^^ only differ in the ordering 
of the expectation values and in that certain moments of 
-j^NPA j^j,g jjQ^ included in x- Still these similarities are 
important to relate results about the NPA method ^^pa 
to the modified moment matrix x- 

PPT Tsirelson bounds & Negativity estimation. — Next 
we demonstrate how the introduced technique helps to 
compute Tsirelson bounds, i.e., extremal quantum val- 
ues of a Bell inequality, or to provide device-independent 
lower bounds on the amount of entanglement. 

The structure of x[p] = xl^i'"] as given by Eq. ([2]) 
provides a device-independent characterization of quan- 
tum correlations: A given probability distribution P can 
only originate from a quantum state if one can find suit- 
able parameters u such that x[p] = xi-Pj""] ^ holds, 
otherwise P is non-quantum. Consequently this super- 
set characterization can be used to derive bounds on the 
optimal quantum value of a Bell inequality, i.e., on a lin- 
ear function of the observed data, generally parsed as 
I ■ P = J2 IabxyP{a, b\x, y), by solving max{/ • P \ x[p] = 
x[^;"] ^ 0}- Increasing the level £ in x[p] generates a 
hierarchy of sharper bounds which, due to the similarity 
between the x[p] matrices and the NPA ones, converges 
to the commutative Tsirelson bound in the limit. Note 
that this optimization is a semidefinite program which 
can be solved efficiently [4^ . 
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In comparison with NPA the advantage of the addi- 
tional structure x = XAb is that one can now incorpo- 
rate further constraints. For instance, one could ask for 
a similar Tsirelson bound if the underlying state is PPT 
by including the additional constraint ^ Oi 



max / • P 

P,u 



(3) 



s.t. 



X[p\ = x[P,n] > 0, x[pf^ > 0, xWtr = 1. 



By this method one again obtains an upper bound to 
the true PPT Tsirelson bound, which converges to the 
related commutative bound in the limit of large levels i, 
see appendix for details. 

Next, let us show how to estimate the negativity [iif . 
defined via the sum of negative eigenvalues Ai of the par- 
tially transposed state as N[pab] — J2\ <o 
the following we employ its variational form which reads 
as iV[pAB]=min{tr[cr_]|pAB= (T+-cr-,cr±^ > 0}. Using 
the properties of the moment matrix, one can readily op- 
timize over a larger set: The constraint p = — cr_ is 
relaxed by x[p\ = x['^+]~x['^-]^ while o'^'* > translates 
to x[(T±]-'"'* > 0. If one observes a certain violation of a 
Bell inequality I ■ P = v, a lower bound on the negativity 
of pab compatible with this observation is given by 



mm X 

P,U,P± ^u± 



s.t. 



x[p] = 



--x[P,u]=xW+]-xW-]>o, x[p]tr 

^x[P±,u±f'^ >0, I-P^v. 



(4) 
= 1, 



At this point we like to stress that both optimization 
problems presented above admit a natural generalization 
to the multipartite scenario using PPT mixtures (which 
include biseparable states) and the genuine negativity as 
a measure for genuine multiparticle entanglement [33|. 
Further details and the explicit semidefinite programs are 
given in the appendix. 

Examples. — Let us start with the Clauser-Horne- 
Shimony-Holt (CHSH) inequality [il], where each party 
has two possible settings x,y G {1,2} yielding binary 
outcomes a, b. Using correlation terms (X^Yy) = P{a — 
b\x, y) — P{a 7^ b\x, y), the inequality reads as 

LHV 

/cHSH = {XiYi) + {X1Y2) + {X2Y1) - {X2Y2) < 2 (5) 



for any local hidden- variable model (LHV), while quan- 
tum mechanics allows a maximum of /chsh = 2\/2. 
Since every separable state fulfills the LHV bound |4J], 
any violation Jchsh > 2 signals entanglement of the un- 
derlying quantum state pab- By solving Eq. Q we can 
now provide a quantitative statement in terms of the min- 
imal negativity that the underlying state pab must pos- 
sess. Specifically, our numerical results lead to the sharp 
bound N[pAB |/CHSH = v] > {v - 2)7(4^2-4); the re- 
sulting plot and a more detailed discussion, also about 
the other examples, can be found in the appendix. Note 



that this recovers the known result that PPT states must 
necessarily satisfy the CHSH inequality (isj . 

As a second example, we consider the Bell inequality 
^3322 < [4^, 1131 where each party can perform three 
possible dichotomic measurements as indicated by the 
subscripts. For < w < 0.25, our numerical solution of 
Eq. Q gives N[pAB 1^3322 =v]>2v. As with the CHSH 
inequality, we observe that the inequality can only be vio- 
lated if the underlying state is NPT entangled [precision: 
10^^]. The other interesting region is close to the max- 
imum possible quantum value. A Bell state can reach a 
violation of 13322 = 0.25 [i^, however the best known 
upper bound is given by I§^22 ^ 0.25088 and there ex- 
ist infinite-dimensional states which can asymptotically 
reach this value [i^. As we see more detailed in the 
first plot of Fig. [U if 13322 > 0.25 the negativity bound 
satisfies N[pab] > 1/2. Because the negativity of any 
(g) state is upper bounded by 1/2 for all D > 2, 
this certifies that such violation is achievable only with 
local Hilbert space dimension d > 3, or, in other words, 
that /3322 < 0.25 serves as a dimension witness [l3| for 
qutrits. 

This insight inspired us to investigate also the very 
first Bell inequalit y u sed as a dimension witness [17| . 



namely, /2233 The corresponding neg- 

ativity bounds are shown in the second plot of Fig. [T] 
At a violation of w = l/\/2 - 1/2 « 0.2071 the neg- 
ativity bound is 1/2, therefore a larger quantum vi- 
olation shows that the underlying state is at least 3- 
dimensional — this also certifies the heuristic qubit bound 
of /2233 presented in Ref. [l3|- Also, for the maximal 
violation at I^^^ = (-\/ll/3 — l)/3, our numerical op- 
timization gives a negativity bound that differs from 
that of the optimal, non-maximally entangled state [Hlj 
ipx cx |00) + X 111) + |22) by less than 5 x 10"^ 

Let us point out that the best PPT Tsirelson upper 
bound that we computed for 12233 does not yet coincide 
with the LHV bound, likewise for a few other Bell in- 
equalities involving a small number of measurement set- 
tings. However, for the majority of known facet-defining 
Bell inequalities involving four dichotomic measurement 
settings per party |52l454| . the respective bounds were 
found to be vanishing within numerical precision. For 
details, see Tab. |T] of the appendix in which we list the 
best (non- vanishing) values for some of these inequalities. 

We also considered examples involving more than two 
parties, where one is typically interested to verify genuine 
multipartite entanglement. This strongest form of mul- 
tiparticle entanglement can be detected from observed 
correlations by violatin g a device- independent entangle- 
ment witness (DIEW) [l^. By a similar method as in 
the bipartite case, we investigated the minimal amount 
of multipartite entanglement (measured by the genuine 
negativity [s^]) needed to violate the DIEWs 132 and 133, 
where each party has respectively two or three dichotomic 
measurements [l^ . [ssj . Let us point out that 132 is the 
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FIG. 1: Negativity bounds for given violations of Bell inequal- 
ities obtained by numerically solving Eq. ([4]) for different lev- 
els of the considered moment matrix. The upper plot shows 
the I3322 inequality [i^ for violations close to the maximum, 
while the lower one corresponds to the 12233 inequality [46| . 



Svetlichny inequality [56[, and thus also demonstrates 
genuine nonlocality. From the bounds we computed (see 
appendix), we can obtain information about the type of 
entanglement resp onsible for given violations, in similar 
spirit to Ref. [14| . For instance, since the genuine neg- 
ativity of any state of the three-qubit W^-class [H^l is 
bounded by \/2/3, one verifies that violations close to 
the maximum of these DIEWs can never be achieved by 
such type of entanglement. Moreover, our bounds show 
that these DIEWs can never be violated by states which 
Using similar arguments as pre- 
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are PPT mixtures 

sented in Ref. 15|, this result can even be extended to 
the n-partite witnesses /„2 and /„3. This implies that, 
apart from a quantification, the generalization of PPT 
Tsirclson bounds to the multipartite case provides a way 
to approximate the set of biseparable correlations in the 
presence of more than three parties [l^l ■ 

Finally, there are also other questions for the multipar- 



tite case. At last we computed the maximal violation of 
the tripartite Bell inequality Is5 < 3 [47| for states which 
are PPT for all bipartitions. We find that it is bounded 
by 3.0187, which shows that the example of Ref. |58|, 
optimally violates the tripartite Peres conjecture via this 
inequality. 

Conclusion. — We have presented a versatile tool to 
characterize quantum correlations in the bi- and mul- 
tipartite case when the underlying state has additional 
properties such as a positive partial transpose. The cru- 
cial observation is that any matrix of moments, which 
originates from a quantum state, can be seen as the re- 
sult of a local map acting on the underlying quantum 
state if the moments are ordered properly. This key idea 
finally allows us to provide device-independent bounds on 
the entanglement present in the measured system. More- 
over, the bounds obtained can provide further informa- 
tion about the underlying state space dimension and the 
type of entanglement involved. Finally, the above obser- 
vation also facilitates the characterization of biseparable 
quantum correlations involving arbitrary number of par- 
ties. 

For future work, we believe that our method can be ex- 
tended to bound, in a device-independent manner, other 
entanglement measures. Clearly, it will also be inter- 
esting to investigate how our technique can be used in 
conjunction with other separability criteria, or applied 
in the closely-related steering scenarios [H^ . 
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Appendix 

Examples of constraints on the matrix of moments 

As an example of how the relations 1) — 3) from 
the main text lead to the general structure of x[p\ 
as given by Eq. we consider the case where the 
map Ka is applied to a single system pA, i.e., x[p\ = 
Sij I*) {i\^^[PAA^jAi]. Suppose that the first setting x 
has three outcomes {0,1,2} whereas the second one, 
here labeled as x', has only two outcomes {0, 1}. Using 
{Ai\ = {l,Mo|^,Mi|^, Afoio:'} one gets 



Mil 



itAfo 



M, 



\x' 



01x^11^ 



Ml 

0\x 



(6) 



PA 
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where pA indicates that we still must take the expectation 
values. Via the listed properties 1) — 3) one obtains 



x[p] 



tr(p) P[Q\x) 

P(0|x) P(0|x) 

P{\\x) 

P{Q\x') Ui ~ iU2 



P{l\x) P{Q\x') 

ui + iu2 

P{l\x) U3 + iU4 

U3 — iw4 P{Q\x') 



(7) 



For the diagonal entries one employs tr[pyiAfo|xAfo|2;] ~ 
ItIpaMq^x] = P(0|x), while the zero entries occur by the 
projection identity MQ|j.Afx|3. = 0. Since the expectation 
value of M^yj-M^yj,, is not directly accessible we can only 
set it equal to a general complex entry ui +iu2 using real 
coefficients u. 

Note that whenever the underlying state p is normal- 
ized, the first term of x[p\ is fixed to be tr(p) = 1. In the 
case where the underlying operator does not need to be 
normalized, however, as it happens for instance with a± 
in the negativity estimation by Eq. ([3]), this entry is not 
fixed a priori. It is thus given by an unknown variable 



The multipartite scenario 

In this section we define more precisely moment ma- 
trices for an arbitrary number of parties, and present 
multipartite optimization problems that compute the 
Tsirelson bounds for PPT mixtures and the device- 
independent quantification of entanglement in terms of 
genuine negativity. 

In a multipartite scenario one can ask, in analogy with 
Eq. ([3]), for a bound on the observed correlations if the 
underlying state is a PPT mixture [33| . A n-partite state 
p is said to be a PPT mixture if it can be written as a con- 
vex combination p ~ '^„^PmPm of states pm which are 
PPT with respect to different bipartitions to C {1, . . . , n} 
of the n subsystems. Since any state that is separable 
with respect to a chosen bipartition m is also PPT ac- 
cording to this splitting, the set of biseparable states is 
included in the set of PPT mixtures. Hence if one verifies 
that a given state is not a PPT mixture, one automat- 
ically certifies that it is not biseparable, and thus, by 
definition, genuine multipartite entangled. 

The definition of moment matrices as given in the main 
text naturally extend to this n-party scenario by apply- 
ing a local CP map to each subsystem s = 1, . . . , ri, 
i.e., p I— x[p\ = ®s=i^s[/o]- By the local structure of 
this transformation one obtains: For any PPT mixture p 
the resulting matrix of moments can be decomposed as 
X[p\ = Y.mPmX[Pm] with x[Pm] > and x[Pm]'^" > 0, 
with m referring to the bipartition on the output spaces. 
Similar properties as given by 1) — 3) in the main text, 
constrain the general structure of x[p\ to be xI-Pj""] f-^d 
one readily obtains a superset approximation for correla- 
tions that can be attained via PPT mixtures. Thus if one 



is interested in the optimal values of a linear expression 
like / • P, where P is generated by a PPT mixture, one 
obtains a bound by solving 

max IP (8) 

S-t., X[P] = X[P,u] = ^x[PmPm] ^^x[Prn,uJ^, 
X[p]tr 1, X[PmPni] > 0, xfPmPm]'^* > Vto. 

Note that in this formulation we included already the 
probabilities pm into the matrix x[Pm,Um\ such that 
x[Pm,Mm]tr = Pni in this case. We show later in the 
appendix that even this multipartite extension converges 
in the limit of an infinite number of moments in x[p]- 

For the multipartite equivalent of Eq. (|4]), let us first 
remind that the genuine negativity [33[ is a computable 
measure of genuine multipartite entanglement which re- 
duces to the negativity in the bipartite case. It is given 

by No[p] = min{p,„,p„} X^mP™^™!/'™] where de- 
notes the negativity with respect to bipartition m and 
the minimization runs over all possible valid state de- 
compositions of the density operator p — J^mPmPm ^■ 
In analogy to Eq. (|4]), if one observes a value of / • P = f , 
a lower bound on the genuine negativity compatible with 
this violation is given by 



(9) 



s-t- X[P] = ^XiPmPrn], I ■ P = V, x[p]tr = 1, 

XiPrnPm] = xW^i] ~ XW^] > 0, X^til^"' > ^TO. 



Further examples and additional information on the 
bounds of negativity/ Bell-inequality violation 

In this part we present more details on the negativity 
bounds mentioned in the main text together with some 
further examples of PPT Tsirelson bounds. 

First note that any explicit quantum state po which 
attains a certain value u of a Bell inequality I ■ P ~ v 
provides an upper bound on the minimal negativity com- 
patible with this violation, i.e., N^in[p\I ■ P = v] = 
minA^[p|J- P = v] < N[po]. Furthermore, since the 
negativity is convex and invariant by adding local auxil- 
iary states, the minimal negativity is a convex function 
in the amount of violation v, that we shall denote by 
f{v) = A^min[p \ I ■ P ~ v]. If one uses the moment matri- 
ces x[p\ with increasing levels £ one obtains lower bounds 



^ The dual of the given measure, i.e., the way how it is de- 
fined in Ref. [H, would be Ng{p) = - min{tr[piy] |Vm : W = 
Pm + Qm™, 1 > Qm > 0, Pm > 0}. Hcnce One sees a missing 
constraint 1 > Pm as compared to its original definition, but 
this does not alter the properties of this measure, i.e., zero for 
biseparable states, full LOCC monotone, convexity and invari- 
ance under local unitaries. 
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fe{v) ^ fiv) with increasing accuracy. However, if one 
finds in the bound fi{v) an interval v S [a, b] such that 
fji{v) is hnear and where the endpoints v ~ a,v ^ b are 
known to be attainable by explicit quantum states then 
fi{v) = f(v) is a tight bound of the minimal negativity 
for V G [a, b]. 

As we see in more detail in Fig. [31 for both, the 
CHSH inequality and the 13322 inequality for a value be- 
tween [0,1/4], the negativity bounds that we obtained 
by solving Eq. (|4]) correspond to straight lines for the 
level i = 3. Specifically, for the CHSH inequality, this 
negativity bound is a straight line joining the points 
(/chsh = 2,N[p] = 0) and {hhsh - 2^/2, Af[p] = 1/2) with 
largest numerical deviation « 7 x 10~^ among all the 
computed instances. Likewise, for the 13322 inequality, 
it is a linear bound connecting the origin and the coor- 
dinate (73322 = 1/4, = 1/2) with largest numerical 
deviation « 3x 10" . Since both endpoints /chsh = 
and /3322 = 1/4 can be achieved with a maximally en- 
tangled two-qubit state, having negativity N[p] = 1/2, 
we thus arrive at the sharp negativity bounds presented 
in the main text. 

Negativity bounds CHSH 




2.83 



Negativity bounds I3; 




0.1 0.15 
Bell vioiation 



FIG. 2: Negativity bounds for CHSH and the J3322 Bell in- 
equality as given by Eq. Q for different levels of the moment 
matrix considered. 



The amount of genuine negativity necessary to achieve 
different violations of the 132 and 133 DIEWs were com- 
puted according to Eq. ©. For biseparable states these 
expressions are bounded by 132 < 4, 133 < 6-\/3, while the 
maximal quantum values are I^'^^ = 4-\/2, 133''^ = 9\/3 
respectively 12, 5^ 56 1. The results for different hier- 
archy levels are shown in Fig. [3l Here we refer to par- 
tial levels such as level 2-1- to denote moment matrices 
X which were constructed with all terms involving 2-fold 
products of local measurement operators, as well as some 
3-fold ones. Similarly to the CHSH and 13322 inequalities, 
the plot for the Svetlichny inequality 132 at level i = 3 
is a straight line up to numerical precision. Since the 
value 732 = 4-\/2 is achievable with a thrce-qubit GHZ 
state with genuine negativity Na[p] = 1/2, we obtain 
the following tight bound for the minimal genuine nega- 
tivity compatible with a Svetlichny inequality violation, 
Ng[pabc\I32 =v]>{v- 4)/(8V2 - 8). 

From Fig. [3] we also note that no violation of these 
DIEW is possible with PPT mixtures; hence the bound 
for biseparable states coincides with the one for PPT 
mixtures. Since any DIEW with this property can be 
found by using Eq. ([5]), and this hierarchy applies to an 
arbitrary number of parties, this means that our tech- 
nique can be used to find DIEWs of this kind for an 
arbitrary number of parties, in contrast to the numerical 
three-party technique presented in Ref. 



12 



Furthermore, as pointed out in the main text, since any 
three qubit state of the W-class [32| satisfies ^ Ng{pw) < 
a/2/3 w 0.471 , our result shows that violations v > 5.563 
and V > 15.36 of the J32 and 133 inequalities can never be 
obtained by any such thrce-qubit states. Therefore one 
gains even further information from the achieved bound 
about the underlying type of entanglement, in similar 
spirit to Ref. 



Finally, note that apart from the bipartite Bell inequal- 
ities mentioned in the main text, we have also computed, 
via Eq. ([3]), upper bounds on PPT Tsirelson bounds 
for 175 facet-defining Bell inequalities involving four di- 
chotomic measurement settin gs p er party. These inequal- 



ities were taken from Refs. |52| - |54| . Interestingly, only 



three of these inequalities give a PPT Tsirelson upper 
bound which differs from the respective LHV bound by 
more than 10 whereas the numerical precision of the 
computation is of the order of 10~^ - 10"^". In Tab. U 
we summarize the best bound that we have computed for 



^ Since Nq is convex, this optimization can be performed over pure 
three-qubit states \'il)w) of the W-class. Note that the genuine 
negativity of a pure state is just the minimum bipartite negativity 
for all bipartitions m. Because Nq is furthermore invariant under 
local basis changes one can employ the Ll/-equivalent standard 
form of a pure three-qubit state of the VF-class [stII , which leaves 
a straightforward optimization. The bound is saturated by the 
VF-state. 
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Genuine negativity bounds Svetlichny I32 




0.5 



0.4 - 



1=1 



4.6 4.8 5 
DIEW violation 

Genuine negativity bounds I33 



not improve further even if one includes higher moments. 
However, in all our considered non-trivial PPT examples 
the respective bounds sharpened if we considered higher 
levels. 

Bell inequality PPT Tsirelson bound I Matrix size 



/chsh < 2 [4^ 

/3322 < [4^ 
/2233 < [46] 

^6 < [52] 
Ji422 < [52] 
4422 < [52] 



2 


< 1.3559 X 10" 

< 7.6754 X 10" 

< 9.7947 X 10" 

< 2.8531 X 10" 



1 3^ X 3^ 

2 10^ X 10^ 
3+ 45^ X 45^ 
2+ 31^ X 31^ 
2+ 29^ X 29^ 
2+ 27^ X 27^ 



/S5 < 3 [43 



< 3.0187 



T X T 



TABLE I: Upper bounds on the maximal possible violation 
by PPT states for different bipartite Bell inequalities via the 
method given by Eq. ((3]). Here matrix size refers to the di- 
mension of the moment matrix x[p] involved in the computa- 
tion whereas I stands for the corresponding level of the mo- 
ment matrix (e.g., level 2-1- means that it contains all 2-fold 
products of local measurement operators, but also some terms 
which only appear when considering 3-fold products of local 
measurement operators). The last inequality corresponds to 
the tripartite case which is PPT for all bipartitions. The pre- 



cision of the computation is of the order of 10 



10" 



6.5 



7 7.5 8 

DiEW violation [in units of VS] 



8.5 



FIG. 3: Genuine negativity as a function of the observed vio- 
lation of the tripartite Svetlichny 732 and DIEW /as inequali- 
ties. Since any three-qubit state of the VK-class has a genuine 
negativity < 0.471, larger violations certify that the underly- 
ing state is not of this SLOCC class. 



these three inequalities as well as for other Bell inequali- 
ties mentioned in the main text. This table also includes 
the respective Tsirelson bound for the tripartite setting 
of 755 with respect to states that are PPT for all biparti- 
tions. This optimization problem is like Eq. ([3]) with the 
tripartite moment matrix and having a PPT constraint 
> for each bipartition. 

Via this numerical investigation we hoped to find a 
counterexample to the bipartite Peres conjecture; a PPT 
state which could violate a Bell inequality. This per- 
spective was triggered by the NPA hierarchy (or using 
also the modified moment matrix) to compute standard 
Tsirelson bounds. Although this method is only guar- 
anteed to be complete in the limit of an infinite num- 
ber of moments, it is important to stress that there are 
many known instances where one could stop the hierar- 
chy already earlier, since one has already reached the true 
Tsirelson bound. This is certified by a special rank prop- 
erty of the solution [4l[ and means that the bound does 



Statement of real states and measurements 

In this section we show that the underlying quan- 
tum state and measurement operators can generally 
be assumed to be real when considering the device- 
independent quantification of entanglement in terms of 
(genuine) negativity, i.e., there exists an equivalent real 
construction having the same (or less) amount of (gen- 
uine) negativity. This extends the result of Rcf. [60[ 
which already proves that probability distributions aris- 
ing from quantum theory in a Bell-type experiment can 
always be reproduced using only real states and real mea- 
surement operators. 

This real property helps in the numerical implemen- 
tation of the bi- or multipartite programs, Eqs. (I3]),(|ll) 
and Eqs. (|8]), (O respectively, since it provides a notable 
parameter reduction in the optimization problems. This 
reduction becomes even greater in the presence of addi- 
tional symmetries. We shortly comment on this param- 
eter reduction at the end of this section. 

Proposition 1. Any n-partite probability distribution 
P having a quantum representation with respect to den- 
sity matrix p and measurement operators Aj, Aj, . . . , A^ 
(with the superscript labeling the party), also has a real 
quantum representation, i.e., a representation in terms of 
a real-valued quantum state p = fF and real-valued mea- 
surement operators A\ = Al'^ , — A^^ , ■ . ■ , A"^ 
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having an equal or lower amount of (genuine) negativ- 
ity, Nm[p] = Nm[p\ for any bipartition m and Ng[p\ < 
Nq[p\. Furthermore one has: 

a) If p is PPT across a bipartition m then p can even be 
assumed to be PPT invariant across m, i.e., p~ p^"^ = 
p^ . Likewise, if p is a PPT mixture then p can even be 
assumed to have a mixture of real PPT invariant states, 

P= Y^PmPm with Pm = P^ = Pm for all TU. 

b) If the observed distribution is invariant under ar- 
bitrary exchange of the parties, P = V(7r)Py(7r)^ for 
all possible permutation tt, then the underlying state and 
measurements can further be assumed to be permuta- 
tionaly invariant, i.e., V{'K)'pV{'iTy = P for all n and 
A\=Al^... = A-. 

Proof. For the first statement we construct for any set 
of projectors , s = 1, . . . ,n and any state p, anotlier 
real state p = fF and real projectors A'^ ~ A^-^ such that 
tr[pA^ (g) (g) . . . (g) A"] = tT[pA^ (g) A'^ g) . . . (g A"] with 
Nm[p\ = Nm[p] for any bipartition m and Ng[p] < Ng[p]. 

Let us start with the projectors. For any chosen basis 
we can decompose the matrix A = A^ -\- iAi into a real, 
symmetric part Aj. = R€(A) = Aj and an imaginary, 
anti-Hermitian part Ai ~ Im(A) = ^Aj. To this matrix 
we now associate the real and symmetric matrix 



1® Ar 



Y ®Ar 



(10) 



where Y = iay. This represents the well-known isomor- 
phism between Hermitian and symmetric matrices, e.g., 
Ref . . Note that the operator A acts on the enlarged 
Hilbert space — 'Ha'®'Ha with 'Ha' — 'C^ ■ Since the 
original A satisfies the projection identity A^ ^ A this 
provides the relations Af — A'f = Ar and ArAi -f AiA,. = 
Ai and henceforth the projection identity for A by 

A^ ^ {1 g) Ar + Y g) Ai){l g)Ar + Yg)Ai) 

= lg){Al- A\) + r ® {ArAi + ArAi) = A. (11) 

We employ this construction for all projectors A'' ^ A'' . 

Next let us define the appropriate extension of the den- 
sity matrix p to p = 'Pa'^AxA'^A'2...a> A^ where the sub- 
scripts label the party to which the Hilbert space belongs. 
However, for ease of reading, we shall use a different or- 
dering of the Hilbert spaces (namely, the primed ones fol- 
lowed by the unprimed ones) PA'^A'^...A',^A^A.i...Ar, = T^{p), 
which could be undone p = V'pA' ...A'^Ai...A„y by an ap- 
propriate permutation V on Hilbert spaces. More pre- 
cisely, this operator is given by 



PA;Ai...A;AiA....A„ - 2-" Re [(1 - iF)®" ® p] 



Let us first show that p is positive semi-definite 
and has trace one as required for any density opera- 
tor. To this end, first note that the right-hand-side 



of Eq. ((12)) can be rewritten as a sum of two orthogo- 
nal parts, namely, 2~"~^(1 — iy)®" (g) p and its trans- 
pose, where the orthogonality follows from the fact that 
(1 — iY){l + iY) = 0. Since the set of eigenvalues re- 
mains unchanged under transposition, this means that 
the eigenvalues of p are precisely one-half of the respec- 
tive eigenvalues of 2~"'(1 — iF)®" (g) p, but with twice 
the multiplicity. However since (1 — iY)/2 is just a rank- 
one projector one has that the non-vanishing eigenval- 
ues of 2~"(1 — zF)®" g) K for any Hermitian matrix 
K are precisely those of K. Thus, we obtain p > 0, 
i|p!|tr = ||p||tr = 1 and ||7?.(p^™)||tr = ||p^'"||tr for any 
bipartition m. 

Next we need to show that these choices indeed pre- 
serve the expectation values. This is best seen using the 
state given by Eq. (|12p and the measurement form of 
Eq. ([To]) together with the following identity 



tvA' [(1 ~iY)g} p {Ig} Ar + Y g) Ai)] 

= tvA' [1] pAr+ tvA' [-iYY] p A, 
+ liA' [Y] pA^+ tiA' [-iY]pAr 
^2pA, (13) 

which follows from Y^ ^ —1 and tr(y) — 0. Apply- 
ing such identities when tracing out each auxiliary space 
yields 

tT[pA^ (g . . . ® 1"] =2"" Re {tr [F(l - iF)®" g) pVx 
ilg)Al + Yg)Al)g)...]} 
=Re (tr[p A^g)...g) A"]) 
= tiip A^ g) ...g)A"]. (14) 

Next, we prove that the negativity remains constant 
Nm[p] ~ Njn[p] for any bipartition m. Thus we need the 
partial transposition of p with respect to m. For simplic- 
ity, we now consider partial transposition with respect 
to A; the general treatment is completely analogous. To 
this end, we remind that the partial transposition is a lin- 
ear operation and thus from Eq. ([T2|) and the fact that 
Y^ = —Y, we have 

iAtXM...A„ = 3^<B [(1 + lY) ® (1 - zF)®"-! ® p] , 
= C/7^[p^■''l]f/^ (15) 

where U = g) g) lAi...An is a unitary matrix and 
we have made use of the fact that axY ax = —Y. Since 
a unitary does not change the eigenvalues and using the 
previously mentioned invariance of the trace-norm we get 
llp^^lltr = \\n{p^^)\\tr = Wp^^Wtr, and hence the state- 
ment that the negativity is unchanged. This statement 
holds for an arbitrary bipartition m. Finally, suppose 
that p = 'YliPmPm is the optimal decomposition for the 
original state p in the definition of the genuine negativity, 
i.e., A'g[p] = '}2iPmNra\prn\- Notc that because the map 
p M- p is linear and preserves positivity one knows that 
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P ~ '^mPmPm IS sl valid dccomposition of p into states 
attributed to different bipartitions m. Since tfie genuine 
negativity is defined by the minimum over all such de- 
compositions, and because Nm[pm\ — Nm[pm\ one read- 
ily obtains Ng[p] < J^PmNmlPm] = N[p], which proves 
the first part of the proposition. 

Henceforth, we shall work with the real state p and the 
corresponding real measurement operators Af = Af^ , 
s = 1, . . . , 71 that can produce the same correlations. As- 
suming this, statement a) is obtained as follows: If the 
original state p is PPT across a bipartition m, then the 
real state p has to be PPT as well. Therefore p^"* > 
can be considered as another quantum state, which, sim- 
ilarly would produce the correct observations since all 
measurements are real, 

P = tT[pAl ® ® . . . ® A^] 

= tr[p(li®l2^...®l^)^"] 

= tiif"^ Al (g) A^j (g) . . . (g) Al] . (16) 

Here, between the first and the second line we trans- 
posed the measurement operator of all parties belonging 
to the bipartition m. Thus also the equal-weight mixture 
p*'™ — {p + p^'")/2 is a real state that is PPT invariant 
with respect to m and that yields the same correlations. 
A similar statement holds for PPT mixtures with partial 
transposition applied to each state pm in the decomposi- 
tion. 

Statement b) follows from a similar argument. For 
simplicity, we provide the proof below only for n = 2 
parties; the generalization to an arbitrary number of par- 
ties is analogous. Let us suppose that p = pab and that 
this state as well as Ai , 13 j are not yet symmetric under 
exchange of the parties.'^ Let us add appropriate local 
auxiliary states |0) , |1) (on system A' , B' respectively) 
to signal whether the state is swapped or not, and con- 
sider the convex combination given by 

pff = I [pAB®\0l)A'B'm + VpABV^®\10)A'B'm] 

= 1[pab®\01)a'B'{01\ 

+V{pab®\01)a'B'{01\)V^ (17) 

where V — V{7r^g) now denotes the swap operator on 
A = A A' and B = BB' . Similarly for the measurements 

such that one finally can check that this is a valid quan- 
tum representation with 

P = [P + l/pyt]/2 = tr[p^Jj^^lf ® If . (19) 



^ At this point we can already assume that the local Hilbert spaces 
are isomorphic Ha — Hb by an appropriate embedding of a 
possibly smaller space in higher dimension. In this way the swap 
operator V = V{iTAg) is well defined. 



From the structure given by Eq. pT)) one observes that 
yO^g is indeed invariant under the swap operator V. Fi- 
nally since the negativity is convex and symmetric (6l| 
one obtains N[p~^] < N[pab], while the LOCC mono- 
tonicity by measuring the primed systems first provides 
A^[p^J^] > N[pab] such that the negativity is indeed in- 
variant. □ 

Let us briefly remark how Proposition [T] reduces the 
number of free parameters for the semidcfinitc programs. 

Firstly, note that since the underlying state and both 
measurements can be chosen to be real, the matrix of mo- 
ments can be assumed to be symmetric x[p] = x[p]'^ ■ For 
the multipartite PPT mixture question of Eqs. ([H]),®, 
this property can be applied to every bipartition, i.e., 
x[Pin] = xiPm]'^ ■ Concerning the negativity finally, this 
property can be further imposed for the operators ap- 
pearing in the variational formulation, since a solution 
given by = cr+ — a~ with {cr^)'^'^ > would imply 
the alternative solution pm = Pm — ("'m)"^ ~ i'^m)'^ with 
ii'^m)'^]'^"' = [('''m)"^'"]^ ^ 0; ^ud hcucc also its equal- 
weight mixture. Therefore in the negativity estimation 
as given by the semidefinite programs of Eqs. (H},®, we 
can set additionally that x[Pm] = x[Pm]"^ and x['^j«] ~ 

The statement a) also simplifies the computation of 
the respective PPT Tsirelson bounds as given by Eq. ([3]) 
or Eq. Since the state can be assumed to be real 

and PPT invariant, this provides the symmetry x[Pm] = 
x[Pm]^ = x[Pm]^" > for all m. 

Finally, whenever one considers a symmetric Bell in- 
equality [ei], i.e., satisfying labxy = hayx in the bipar- 
tite case, one can impose this symmetry also for the cor- 
responding distribution P{a,b\x^y) = P{h,a\y,x), such 
that b) of Proposition [1] gets relevant. If both local 
maps in the construction of x[p] have an equal number 
of moments, then x[p] can be assumed to be invariant 

Vabx[pWab ~ -^['^l under the swap operator Vab- This 
symmetry in particular helps in the higher levels of the 
hierarchy. 

Convergence of respective PPT characterizations 

As emphasized in the main text, the method to approx- 
imate PPT Tsirelson bounds via Eq. ([3]) converges to a 
description of the commutative set in the limit of an in- 
finite number of moments. A similar statement holds for 
the multipartite case and PPT mixtures. To show this, 
we start this section by formulating the commutative ver- 
sion of the respective quantum representation. We then 
proceed to show the equivalence between this commu- 
tative and the tensor product version in the case of a 
finite-dimensional quantum representation and, at last, 
prove the convergence. We stress that this convergence 
follows from the convergence of the NPA hierarchy [4lj . 
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At first let us motivate this distinction between the 
tensor product and commutative version: In quantum 
information we are commonly used to employ tensor 
products for measurements on different subsystems, i.e., 
tr[pA_B^ ^ B]. In contrast, in algebraic field theory for 
instance, measurements on different parts are described 
by commuting operators [^,-6] = 0, both acting already 
on the bipartite state p (therefore we use the tilde) such 
that expectation values become tr[joAi3]. If one identifies 
A = A(^l, B = l^B and p = pab then one sees that aU 
expectation values using tensor products can be recov- 
ered with commuting observables. The converse ques- 
tion, however, is still open and is known as Tsirelson's 
problem 

For clarity, we start with the translation of the bi- 
partite PPT constraint p^^ > in commutative terms. 
We employ the simplified notions Ai,Bj^ . . . introduced 
in the main text for measurement operators, and P = 
tr[. . . Ai . . .Bj . . ] to indicate the functional form of the 
observed probability distribution. 

Definition 1. i) A bipartite probability distribution P 
has a tensor product PPT quantum representation if 

there exists a real state pAB, local real measurements op- 
erators Ai,Bj such that P = tT[pABAi (X) Bj] and that 
PAB = p"ab PPT invariant. 

ii) A bipartite distribution P has a commutative PPT 
quantum representation if there exists a real state p, com- 
muting real measurement operators Ai,Bj, [Ai,Bj] = 0, 
such that P = trlpAiBj] and for all possible products, 
indexed by i ~ (ii, 12, ... , in), A{i) = Ai-^Ai^ . . . Ai^, and 
similarly for B , it holds 

tT[pA{i)B{j)] - tr[pi(*^)B(j)] (20) 

with = (i„,i„-i, . . . ,ii). 

Proposition 2. Any bipartite probability distribution 
which has a tensor product PPT quantum representa- 
tion also has a commutative PPT quantum representa- 
tion. The converse holds at least if the commutative PPT 
quantum representation is finite- dimensional. 

Proof. Let us start with the direction i) to ii). The for- 
mulation i) is already simplified using Proposition [1] If 
we now interpret A = A^l, B = 1®B and p — pab then 
these operators naturally fulfil the requirements of being 
real, that their commutator vanishes and that they pro- 
vide the correct expectation values. If one employs the 
PPT invariance pab = Pa^b ^^'^ ~ '^^^ obtains 
the condition given by Eq. ([20)) via 

tT[pA{i)B{j)] ^ tT[pABA{i)®B{j)] ^ tT[p'^j,A{i)®B{j)] 
= tT[pABiA^,A,^ . . . A,J'^ (E}B{j)] 

= tr[/^AB(^^„A^„_l• . -A,,)® B{j)] 
= tT[pABA{i'^)®B(])] 

= tr[pA{(^)B{j)]. (21) 



This finishes the proof in the direction i) to ii) . 

Let us now turn to the converse. Thus suppose that 
p and Ai, Bj, all acting on the finite-dimensional Hilbert 
space a, are given as in ii). Prom these operators we now 
construct local measurement operators and a bipartite 
quantum state which is PPT; the statement that they 
even have the extra symmetries follows from point a) 
of Proposition [1] This construction is very analogous to 
the corresponding one of the standard Tsirelson problem. 
However, in order to show that the conditions given by 
Eq. (PO)) finally prove that the constructed bipartite state 
is PPT one has to keep track which operators can be 
built up by products of Ai,Bj and linear combinations 
of them, or, more precisely, by the operators set 

Q = {g e ]B(H)|Q = ^c,,i(^)^(.7),cy■ G (D}. (22) 

From the given finite-dimensional commuting Hcrmitian 
operators Ai,Bj one can infer the following decomposi- 
tions 

ii = ®k'HAB,k ® ICk (23) 
= ®k{®iHAM®'^B,ki)®^k. (24) 

and representations 

Q = ®kHHAB.k)®l, (25) 
A^ = ®kAi,k ® 1 = ®k{®iKki «) 1) «) 1, (26) 
Bj = ®kBj,k (®l = ®fe(®a ® B],ki) ® 1- (27) 

These statements can be inferred for instance from The- 
orem A. 7 of Ref. [sH: The original structure given by 
Eq. comes from the decomposition of (the finite- 
dimensional C*-algebra) Q. The tensor product struc- 
ture of HAB.k, included in Eq. ([M)) . originates from the 
same argument as in the standard Tsirelson problem, ap- 
plied to the operators Ai ]^,Bj ]^ acting on 'HAB,k which 
still commute. (In addition to Theorem A. 7 of Ref. [sil 
one needs to know that Bj k are elements of the commu- 
tant of C*-algebra generated by {Ai^k})- 

Since it will become important shortly, let us stress the 
meaning of Eq. (f25|) : It states that any operator of this 
form is in the set Q. Note that this structure includes, 
in particular, all operators (9k[(SiQki] 1 with Qki ar- 
bitrary. Thus, also the very special operators where all 
of these operators Qki are zero except for one particu- 
lar pair, maybe k = k',l = I'. Hence by appropriately 
chosen coefficients Cij[Qk'i'] it is possible to build up 

J2 [Qk'i']Aii)Bij) - Qk'i' ® Ik,, (28) 

with Qk'i' being an arbitrary operator acting on Hab.m- 
Next let us construct the bipartite state. In the follow- 
ing we denote by H^; the projections onto T-Labm ® l^k- 
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From the structure of the Hilbert space and the measure- 
ments one obtains the identity 

tv[pA,Bj]^^tvABM[^^Kk {nkipnM)Ai^M®i3j^ki\- (29) 

kl 

If one defines pmPabm = tr^^ (rifeipllfe/) one can inter- 
pret Eq. as a classical mixture of bipartite states on 
different subsystems. Via appropriate local auxiliary sys- 
tem one can combine this into a single bipartite system 

PAB = y^^PkiPABM <^ \kl, kl) j^,g, {kl, kl\ , (30) 
A. = ^i,,fc/®|A:0^, (fc/|, (31) 
^^■=II4>i®|fcOi3'(fc^h (32) 



such that tr[pA^Bj] = tr[pAsAi (8) Bj] following Eq. (f29| . 
That PAB is indeed a valid quantum state, while Ai,Bj 
are correct operators describing measurements follows 
from their structure. 

Thus we are left to show that pab is PPT, which, by its 
form given by Eq. pop . is equivalent to pAB,ki being PPT 
for all k, I. A finite-dimensional operator p^-^ fc/ > is 
positive semidefinite if and only if for all operators Qm € 
B{Uabm) it holds that tT[p'^% ,^iQkiQli] > [H. In 
the remaining we show that this inequality holds due to 
tT[pQQ^ > for all operators Q G Q and the extra condi- 
tions given by Eq. pO)) . Here it is important to stress that 
Q indeed includes aU operators Qki £ ^(T-Labm) as re- 
marked previously. Using Q G Q with Cij = Cij[Qk'v\ as 
given in Eq. ^ for fixed k' , V , Aki{i) = Ai^^^i ■ ■ ■ A^m 
and a similar shorthand for Bki{j), one obtains 

tr[pOQt] ^Y.c,,c:,tr[pA{t)B{j){A{u)B{v))^] 

= Y.^^^<vtr[pAi^)Aiu^)B{J)Bivy] 

P ^ c.,cL tT[pAiu)Ai^^)mB{vy] 

= X! ^iJ^uvPkl 

X ti-AB.ki[PABMi^kiii)Aki{u)^)'^ «) Bkiij)i3kiivy] 
X iYAB,kAp5iBM^ki{i)Aki(u)'^ ® Bki{j)Bkiiv)'^] 

= Pk'l' trAB,k'l'[PAB.k'l'Qk'l'Qk'l'] ^0. (33) 

Here we employ the following: Because Ai commutes 
with every Bj this property naturally extends to products 
[A{i), B{j)] = 0. Moreover, since the measurement oper- 
ators Ai can be taken real and symmetric, it follows that 
i(i)t = A{i^) and similarly A{k)A{f) = [A{i)A{kyf. 
As indicated the identity given by Eq. ([20]) is employed 
in the third line. Afterwards one uses once more the 
block decomposition as already shown previously, and 



the identity tr[Xy^-^] = tr[X^^y]. The final inequality 
stems from the positivity of the original density operator 
p. This finishes the proof of Proposition [2] □ 

Let us now see how the property of a PPT mixture can 
be cast into a commutative version. The equivalence of 
both formulations in the finite case is shown similarly as 
in the bipartite case. 

Definition 2. i) An n-partite probability distribution 
P has a tensor product PPT mixture quantum rep- 
resentation ij there exists a real n-partite state p, lo- 
cal real measurement operators A\,A^^,..., such that 
P = tr[pAl (g) ^2 ^ . . . (8) A^] and p = Y^mP^nPrn be- 
ing a convex combination of real, PPT invariant states 
Pm = Pm™ '^^^ hipartitions m. 

ii) An n-partite probability distribution P has a com- 
mutative PPT mixture quantum representation if there 
exists a real state p, commuting real measurement op- 



erators A} , Aj 



,Al such that P = \.x\pA\A\ 



and p = PmPm is a convex combination of real 
'PPT states" pm for all possible formal bipartitions 
m, i.e., each state pm must satisfy that for all possi- 
ble products of operators from partition m, indexed by 
j(m) ^ {i\^,il^^...,i'^rj yjiiii ks € m for s = l,...,r, 

A(i('"^) = A'i^A'^I^ . . . A'l^ and similar for products for its 
complement rUc it holds 

tr[p„i(z(™))i(/™^))] =tr[p™i(i(")^)i(/™=))] (34) 



with = (iJ: 



I. ./c,,_ 



Proposition 3. Any probability distribution which has a 
tensor product PPT mixture quantum representation also 
has a commutative PPT mixture quantum representation. 
The converse holds at least if the commutative PPT mix- 
ture quantum representation is finite- dimensional. 

Proof. The direction from i) to ii) holds by similar argu- 
ments as given in the bipartite case. 

For the converse direction one proceeds analogously as 
in the proof of Proposition 2. From the given measure- 
ment operators one can obtain, similar to Eas. (P5)) -(P7 1) . 
a decomposition of the underlying Hilbert space given by 



n = (Ski^iUA.M^^A^.ki 



^nA^M)®^k, (35) 



its measurement operators and respective operator set Q, 
i.e., the set of all operators generated by linear combina- 
tions of products of the measurements 3l| . This shows 
that all measurements A^ from system s act only non- 
trivial on the part T-La^M^ in analogy to the bipartite 
ease as given by Eqs. ((26|) . p7l) . Via the projector H^; 
onto the subspace 'HAi...A„,ki 8) ICk one arrives again at 
the tensor product form 



tiipAl...A]] 

= Y^Pkl trAi...A„.kl[^klAl f 



(36) 



kl 
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with, pkiaki = Pki^Ai...A„,ki = tr^jrEfcipIIfc;]. Via aiixil- 
iary states this classical mixture can be turned into a sin- 
gle multipartite state, i.e., p = J2kiPki^kiv (81 \kl) 
and measurements Af — ki ® 1^0 (''^^l fo'" each sys- 

tem s. 

At last one needs to show that p is a PPT mixture, 
which amounts to verifying that each dki is a PPT mix- 
ture. This decomposition is given by the original expan- 
sion p = ^^7X1 Pm into real (in this case un-normalized) 
states pm which satisfy the property given by Eq. 
If we set prn,ki ~ ^T^Kki^kiPm^ki] then one obtains a de- 
composition 

PkiS-ki = ^^Pm,kh (37) 

rn 

which verifies that aki is a PPT mixture if one can 
show that pm,ki is PPT with respect to bipartition m. 
Since we can interpret Pm,ki as a bipartite system on 
m and its complement to this boils down to the bi- 
partite argument: The condition pj^\i > is equiv- 
alent to ti[p^"^,,iQkiQli\ > for all operators Qki £ 
IB('HAi...A„,fei)- Since any operator Qki can again be writ- 
ten as a linear combination of products of the measure- 
ment operators, the identities given by Eq. p4p verify 

^^iPmM^kiQli] = tj:[pm,kiQkiQli] > 0, where the in- 
equality results from the positivity of pm ■ This concludes 
the proof. □ 

Finally let us state the convergence result. As already 
pointed out this is a corollary of the original convergence 
proof for the NPA hierarchy [U . In this part we use the 
notation x'"^^ to refer to the matrix of moments of level i, 
i.e., having for each party a product of up to £ local oper- 
ators in the construction. We first prove the convergence 
of the PPT hierarchy for the bipartite case. This proof 
then directly generalizes to a PPT hierarchy for a single 
bipartition m in the multipartite case, which finally helps 
to prove the convergence of the PPT mixture in the end. 
Note that the condition x*-^-* [-f,'"]tr = 1 does not appear 
since we explicitly talk about probability distributions. 

Proposition 4. // a bipartite distribution P satisfies 
X^^^[P,u\ x^'^'[-P,^t]^-* = x'-'^HP^uf > via suitable 
u for all levels £, then it has a (potentially infinite- 
dimensional) commutative PPT quantum representation. 

Proof. In the following we use the label s,t,v,w to de- 
note sequences of measurement operators, e.g., s = 
(si, S2, . . . , Sn) with A{s) — Asj^As2 ■ ■ . As^ . However, if 
we employ i,j it refers to a single index, and i = or 
j = should correspond to the identity operator. The 
structure of x^^\ as given by Eq. ([T]) for ^ = 1, is such 
that the matrix elements correspond to expectation val- 
ues 

{st\x^^^\vw) = tr{p[A{s) B{t)] [A{v) ® B{w)]^. (38) 



The identities given by the described partial informa- 
tion can be parsed as {st\x'^^^vw) = {s't'\x'"^^v'w') if 
Ais)A{vy ® B{t)B{wy = A{s')A{v'y B{t')B{wy by 
the properties satisfied by operators, which are listed 
as 1) — 3) in the main text. The two extra require- 
ments provide the additional relations {st\x^^\vw) = 
{vt\x'^^^\sw) = {vw\x^^^\st). 

Next let us summarize the convergence idea of the NPA 
hierarchy [4l|. Suppose that x^^\P,u] = x^^\P,uV-^ = 
X^^^lP^uy > exists for each level i, then there is also 
a limit x^°°^ > which satisfies all the listed identities. 
For this infinite-dimensional matrix one can associate a 
set of vectors {|est)} with {st\x''°°^vw) = {est\evw)- Via 
this set one now defines the measurements and the corre- 
sponding state. The measurements are A, = proj({|est) : 
si = i}) where proj stands for the projector onto the cor- 
responding set. Similarly one defines 13 j = proj ({| est) : 
ti ~ j})i while the state is given by yO = |eoo) (eoo|. As 
shown in detail in Ref. [4lj these measurements satisfy 
Ai \est) = \est) with s = (i,s) and similar for Bj due 
to the linear constraints. These choices reproduce the 
observed expectation values since 

Pij = (eoo|eij) = {eoQ\A,Bj\eoQ) 

= tr[|eoo) (eool A,B,] = tv[pA,B,]. (39) 

Moreover, besides being true projectors, the operators 
Ai.Bj fulfil the projector identity and commute. 

Now let us add the remaining properties. First, let us 
point out that since x'"°°'' = x'"°°''"^ the set {|est)} can be 
chosen to be a set of real vectors. Therefore, the state and 
also the measurements are real symmetric, which implies 
in particular A{s)^ — A{s'^). Thus wc are left to show 
the property given by Eq. (pO|) . However, this follows 
from the extra requirement {st\x'^•°°^vw) = {vt\x^°°^sw) 
since 

tT[pA{s)B{t)] - {eoo\A{s)B{t)\eoo) 

= (eoolest) = (csoleot) = {eoo\ A{s)^ B{t)\eoo) 

= (eoo|i(s^)^(t)|eoo) = tr[pAis^)B{t)]. (40) 

This concludes the convergence proof. □ 

Let us remark that the proof works equivalently in 
the multipartite case if one is only interested in a sin- 
gle bipartition to: If an n-partitc distribution P sat- 
isfies x'-^HP^u] = X^'^^P^uf"' = X^^HP.uf > via 
suitable u for all levels £, then it has a (potentially 
infinite-dimensional) commutative PPT quantum repre- 
sentation with respect to partition to. As for the bi- 
partite case one uses the convergence statement in the 
multipartite case to show existence of a respective state 
p and commuting measurements Aj, . . . , A^. The addi- 
tional linear requirements given by Eq. p4p for to are 
shown similarly as in Eq. (|40l) by using the identities 
xW[P,u] = x^'\P,uV^ = X^'\PM^- 
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Before we continue we like to clarify another point: 
As seen in the proof of Proposition |4] one always ob- 
tains a pure state p = IV') (V'l in the commutative ver- 
sion. However, since any pure PPT state of a tensor 
product Hilbert space is separable and thus possesses a 
LHV model, one could ask whether this implies that any 
(finite-dimensional) PPT state must necessarily satisfy 
all Bell inequalities. This would prove the Peres conjec- 
ture, at least, in the finite-dimensional case. However 
this is not true. Indeed, even if the state in the commu- 
tative version is pure, the related bipartite, tensor prod- 
uct state following the procedure outlined in the proof of 
Proposition [2] does not have to be. More precisely, this 
bipartite state pab is given by Eq. pO|) . where the partial 
trace over JCu in pabm = ^"^Kki^ki IV") (V'l ^ki) results in 
a mixed PPT state. This extra system ICk can be consid- 
ered as a purifying system of each bipartite state pAB,ki 
to Hki 14'), which would then be pure if one translates it 
into commutative terms. Finally note that the purifying 
systems K-k would vanish if the commuting measurements 
Ai, Bj would generate 13 (H). 

At last, we prove the convergence of the characteriza- 
tion of an n-partite distribution P with an underlying 
PPT mixture, which shows that the hierarchy given by 
Eq. ([5]) is complete. 

Proposition 5. // an n-partite distribution P satis- 

fies X^^HP,u] = EX^^HPm,Urn] With X^'^Pm^Um] = 
X^^HPm,Wm]'^* =X^^HPm,Um]'^ > via Suitable Pm,Um 

for all levels £, then it has a (potentially infinite- 
dimensional) commutative PPT mixture quantum repre- 
sentation. 

Proof. From the remark after Proposition|4]we know that 
the conditions on x^^''[^m,Wm] for all i show that there 
exists an (un-normalizcd) distribution Pm which has a 
commutative PPT quantum representation for the bipar- 
tition m. Let us refer to this representation as the state 
am and measurements A^'™, . . . , A^'™- Note that 

here, because the probabilities P^ are not normalized 
but we rather have J2m Pm ~ P, the state pm satisfies 

tr(p,„) = x[P,„,Um]tr- 

From this one can construct a commutative PPT mix- 
ture representation for the given n-partite distribution P 
by using % = ®mHm, P = ®ra<^ra and Al = ®rnAf''^' for 
each system s. This construction directly ensures positiv- 
ity and that the measurements commute. A commutative 
PPT mixture decomposition of the state is then given by 
p = Pm with (un-normalized) states pm = HmpHm 
and Hm denoting the projector onto That these 

states Pm indeed satisfy the relations given by Eq. (p4)) 
follows from the corresponding identities on tT,„ since 
PmAl — (Jm^i for each system s and similarly for prod- 
ucts of more operators. This concludes the proof. □ 



